Moderate mutation rate in the SARS coronavirus genome and its implications

Background: The outbreak of severe acute respiratory syndrome (SARS) caused a severe global epidemic in 2003 which led to hundreds of deaths and many thousands of hospitalizations. The virus causing SARS was identified as a novel coronavirus (SARS-CoV) an

Estimation of male-to-female ratios of mutation rate in primates and rodents

Haldane (1935) developed a method for estimating the male-to-female ratio of mutation rate ($\alpha$) by using sex-linked recessive genetic disease, but in six different studies using hemophilia A data the estimates of $\alpha$ varied from 1.2 to 29.3. Direct genomic sequencing is a better approach, but it is laborious and not readily applicable to non-human organisms. To study the sex ratios of mutation rate in various mammals, I used an indirect method proposed by Miyata et al. (1987). This method takes advantage of the fact that different chromosomes segregate differently between males and females, and uses the ratios of mutation rate in sequences on different chromosomes to estimate the male-to-female ratio of mutation rate. I sequenced the last intron of ZFX and ZFY genes in 6 species of primates and 2 species of rodents; I also sequenced the partial genomic sequence of the Ube1x and Ube1y genes of mice and rats. The purposes of my study in addition to estimation of $\alpha$'s in different mammalian species, are to test the hypothesis that most mutations are replication dependent and to examine the generation-time effect on $\alpha$. The $\alpha$ value estimated from the ZFX and ZFY introns of the six primate specise is ${\sim}$6. This estimate is the same as an earlier estimate using only 4 species of primates, but the 95% confidence interval has been reduced from (2, 84) to (2, 33). The estimate of $\alpha$ in the rodents obtained from Zfx and Zfy introns is ${\sim}$1.9, and that deriving from Ube1x and Ube1y introns is ${\sim}$2. Both estimates have a 95% confidence interval from 1 to 3. These two estimates are very close to each other, but are only one-third of that of the primates, suggesting a generation-time effect on $\alpha$. An $\alpha$ of 6 in primates and 2 in rodents are close to the estimates of the male-to-female ratio of the number of germ-cell divisions per generation in humans and mice, which are 6 and 2, respectively, assuming the generation time in humans is 20 years and that in mice is 5 months. These findings suggest that errors during germ-cell DNA replication are the primary source of mutation and that $\alpha$ decreases with decreasing length of generation time. ^

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

The mutation rate and cancer

The selection of advantageous mutations underlies tumorigenesis. The growth of a tumor is therefore a form of evolution at the somatic level, in which the population is comprised of individual cells within the tumor. Models of tumorigenesis have considered the relative importance of mutation and selection. We show that selection is more important than an increased mutation rate in the growth of a tumor. Some cancers may acquire a “mutator phenotype,” probably leading to faster growth, but mutator phenotypes are not necessary for carcinogenesis.

Mutation rate varies among alleles at a microsatellite locus: Phylogenetic evidence

The understanding of the mutational mechanism that generates high levels of variation at microsatellite loci lags far behind the application of these genetic markers. A phylogenetic approach was developed to study the pattern and rate of mutations at a dinucleotide microsatellite locus tightly linked to HLA-DQB1 (DQCAR). A random Japanese population (n = 129) and a collection of multiethnic samples (n = 941) were typed at the DQB1 and DQCAR loci. The phylogeny of DQB1 alleles was then reconstructed and DQCAR alleles were superimposed onto the phylogeny. This approach allowed us to group DQCAR alleles that share a common ancestor. The results indicated that the DQCAR mutation rate varies drastically among alleles within this single microsatellite locus. Some DQCAR alleles never mutated during a long period of evolutionary time. Sequencing of representative DQCAR alleles showed that these alleles lost their ability to mutate because of nucleotide substitutions that shorten the length of uninterrupted CA repeat arrays; in contrast, all mutating alleles had relatively longer perfect CA repeat sequences.

The high spontaneous mutation rate: Is it a health risk?*

The human mutation rate for base substitutions is much higher in males than in females and increases with paternal age. This effect is mainly, if not entirely, due to the large number of cell divisions in the male germ line. The mutation-rate increase is considerably greater than expected if the mutation rate were simply proportional to the number of cell divisions. In contrast, those mutations that are small deletions or rearrangements do not show the paternal age effect. The observed increase with the age of the father in the incidence of children with different dominant mutations is variable, presumably the result of different mixtures of base substitutions and deletions. In Drosophila, the rate of mutations causing minor deleterious effects is estimated to be about one new mutation per zygote. Because of a larger number of genes and a much larger amount of DNA, the human rate is presumably higher. Recently, the Drosophila data have been reanalyzed and the mutation-rate estimate questioned, but I believe that the totality of evidence supports the original conclusion. The most reasonable way in which a species can cope with a high mutation rate is by quasi-truncation selection, whereby a number of mutant genes are eliminated by one “genetic death.”

Multilocus DNA fingerprinting reveals high rate of heritable genetic mutation in herring gulls nesting in an industrialized urban site.

Genotoxins, such as polycyclic aromatic compounds, are ubiquitous in urban and industrial environments. Our understanding of the role that these chemicals play in generating DNA sequence mutations is predominantly derived from laboratory studies with specific genotoxins or extracts of contaminants from environmental media. Most assays are not indicative of the germinal effects of exposure in situ to complex mixtures of common environmental mutagens. Using multilocus DNA fingerprinting, we found the mutation rate in herring gulls inhabiting a heavily industrialized urban harbor (Hamilton Harbour, Ontario) to be more than twice as high as three rural sites: Kent Island, Bay of Fundy; Chantry Island, Lake Huron; and Presqu'ile Provincial Park in Lake Ontario. Overall we found a mutation rate of 0.017 +/- 0.004 per offspring band in Hamilton, 0.006 +/- 0.002 at Kent Island, 0.002 +/- 0.002 from Chantry Island, and 0.004 +/- 0.002 from Presqu'ile Provincial Park. The mutation rate from the rural sites (pooled) was significantly lower than the rate observed in Hamilton Harbour (Fisher's exact test, two-tailed; P = 0.0006). These minisatellite DNA mutations may be important biomarkers for heritable genetic changes resulting from in situ exposure to environmental genotoxins in a free-living vertebrate species.

Mutation frequency of non-ESBL phenotype SENTRY (Asia-Pacific) isolates of Klebsiella pneumoniae conversion to an ESBL positive phenotype

Extended spectrum β-lactamases or ESBLs, which are derived from non-ESBL precursors by point mutation of β-lactamase genes (bla), are spreading rapidly all over the world and have caused considerable problems in the treatment of infections caused by bacteria which harbour them. The mechanism of this resistance is not fully understood and a better understanding of these mechanisms might significantly impact on choosing proper diagnostic and treatment strategies. Previous work on SHV β-lactamase gene, blaSHV, has shown that only Klebsiella pneumoniae strains which contain plasmid-borne blaSHV are able to mutate to phenotypically ESBL-positive strains and there was also evidence of an increase in blaSHV copy number. Therefore, it was hypothesised that although specific point mutation is essential for acquisition of ESBL activity, it is not yet enough, and blaSHV copy number amplification is also essential for an ESBL-positive phenotype, with homologous recombination being the likely mechanism of blaSHV copy number expansion. In this study, we investigated the mutation rate of non-ESBL expressing K. pneumoniae isolates to an ESBL-positive status by using the MSS-maximum likelihood method. Our data showed that blaSHV mutation rate of a non-ESBL expressing isolate is lower than the mutation rate of the other single base changes on the chromosome, even with a plasmid-borne blaSHV gene. On the other hand, mutation rate from a low MIC ESBL-positive (≤ 8 µg/mL for cefotaxime) to high MIC ESBL-positive (≥16 µg/mL for cefotaxime) is very high. This is because only gene copy number increase is needed which is probably mediated by homologous recombination that typically takes place at a much higher frequencies than point mutations. Using a subinhibitory concentration of novobiocin, as a homologous recombination inhibitor, revealed that this is the case.

Time dependency of molecular rate estimates and systematic overestimation of recent divergence times

Studies of molecular evolutionary rates have yielded a wide range of rate estimates for various genes and taxa. Recent studies based on population-level and pedigree data have produced remarkably high estimates of mutation rate, which strongly contrast with substitution rates inferred in phylogenetic (species-level) studies. Using Bayesian analysis with a relaxed-clock model, we estimated rates for three groups of mitochondrial data: avian protein-coding genes, primate protein-coding genes, and primate d-loop sequences. In all three cases, we found a measurable transition between the high, short-term (<1–2 Myr) mutation rate and the low, long-term substitution rate. The relationship between the age of the calibration and the rate of change can be described by a vertically translated exponential decay curve, which may be used for correcting molecular date estimates. The phylogenetic substitution rates in mitochondria are approximately 0.5% per million years for avian protein-coding sequences and 1.5% per million years for primate protein-coding and d-loop sequences. Further analyses showed that purifying selection offers the most convincing explanation for the observed relationship between the estimated rate and the depth of the calibration. We rule out the possibility that it is a spurious result arising from sequence errors, and find it unlikely that the apparent decline in rates over time is caused by mutational saturation. Using a rate curve estimated from the d-loop data, several dates for last common ancestors were calculated: modern humans and Neandertals (354 ka; 222–705 ka), Neandertals (108 ka; 70–156 ka), and modern humans (76 ka; 47–110 ka). If the rate curve for a particular taxonomic group can be accurately estimated, it can be a useful tool for correcting divergence date estimates by taking the rate decay into account. Our results show that it is invalid to extrapolate molecular rates of change across different evolutionary timescales, which has important consequences for studies of populations, domestication, conservation genetics, and human evolution.

Evidence for time dependency of molecular rate estimates

Long-term changes in the genetic composition of a population occur by the fixation of new mutations, a process known as substitution. The rate at which mutations arise in a population and the rate at which they are fixed are expected to be equal under neutral conditions (Kimura, 1968). Between the appearance of a new mutation and its eventual fate of fixation or loss, there will be a period in which it exists as a transient polymorphism in the population (Kimura and Ohta, 1971). If the majority of mutations are deleterious (and nonlethal), the fixation probabilities of these transient polymorphisms are reduced and the mutation rate will exceed the substitution rate (Kimura, 1983). Consequently, different apparent rates may be observed on different time scales of the molecular evolutionary process (Penny, 2005; Penny and Holmes, 2001). The substitution rate of the mitochondrial protein-coding genes of birds and mammals has been traditionally recognized to be about 0.01 substitutions/site/million years (Myr) (Brown et al., 1979; Ho, 2007; Irwin et al., 1991; Shields and Wilson, 1987), with the noncoding D-loop evolving several times more quickly (e.g., Pesole et al., 1992; Quinn, 1992). Over the past decade, there has been mounting evidence that instantaneous mutation rates substantially exceed substitution rates, in a range of organisms (e.g., Denver et al., 2000; Howell et al., 2003; Lambert et al., 2002; Mao et al., 2006; Mumm et al., 1997; Parsons et al., 1997; Santos et al., 2005). The immediate reaction to the first of these findings was that the polymorphisms generated by the elevated mutation rate are short-lived, perhaps extending back only a few hundred years (Gibbons, 1998; Macaulay et al., 1997). That is, purifying selection was thought to remove these polymorphisms very rapidly.

Mutations as molecular tools : The metabolic-rate dependent molecular clock and DNA barcoding of allied species

Mutation and recombination are the fundamental processes leading to genetic variation in natural populations. This variation forms the raw material for evolution through natural selection and drift. Therefore, studying mutation rates may reveal information about evolutionary histories as well as phylogenetic interrelationships of organisms. In this thesis two molecular tools, DNA barcoding and the molecular clock were examined. In the first part, the efficiency of mutations to delineate closely related species was tested and the implications for conservation practices were assessed. The second part investigated the proposition that a constant mutation rate exists within invertebrates, in form of a metabolic-rate dependent molecular clock, which can be applied to accurately date speciation events. DNA barcoding aspires to be an efficient technique to not only distinguish between species but also reveal population-level variation solely relying on mutations found on a short stretch of a single gene. In this thesis barcoding was applied to discriminate between Hylochares populations from Russian Karelia and new Hylochares findings from the greater Helsinki region in Finland. Although barcoding failed to delineate the two reproductively isolated groups, their distinct morphological features and differing life-history traits led to their classification as two closely related, although separate species. The lack of genetic differentiation appears to be due to a recent divergence event not yet reflected in the beetles molecular make-up. Thus, the Russian Hylochares was described as a new species. The Finnish species, previously considered as locally extinct, was recognized as endangered. Even if, due to their identical genetic make-up, the populations had been regarded as conspecific, conservation strategies based on prior knowledge from Russia would not have guaranteed the survival of the Finnish beetle. Therefore, new conservation actions based on detailed studies of the biology and life-history of the Finnish Hylochares were conducted to protect this endemic rarity in Finland. The idea behind the strict molecular clock is that mutation rates are constant over evolutionary time and may thus be used to infer species divergence dates. However, one of the most recent theories argues that a strict clock does not tick per unit of time but that it has a constant substitution rate per unit of mass-specific metabolic energy. Therefore, according to this hypothesis, molecular clocks have to be recalibrated taking body size and temperature into account. This thesis tested the temperature effect on mutation rates in equally sized invertebrates. For the first dataset (family Eucnemidae, Coleoptera) the phylogenetic interrelationships and evolutionary history of the genus Arrhipis had to be inferred before the influence of temperature on substitution rates could be studied. Further, a second, larger invertebrate dataset (family Syrphidae, Diptera) was employed. Several methodological approaches, a number of genes and multiple molecular clock models revealed that there was no consistent relationship between temperature and mutation rate for the taxa under study. Thus, the body size effect, observed in vertebrates but controversial for invertebrates, rather than temperature may be the underlying driving force behind the metabolic-rate dependent molecular clock. Therefore, the metabolic-rate dependent molecular clock does not hold for the here studied invertebrate groups. This thesis emphasizes that molecular techniques relying on mutation rates have to be applied with caution. Whereas they may work satisfactorily under certain conditions for specific taxa, they may fail for others. The molecular clock as well as DNA barcoding should incorporate all the information and data available to obtain comprehensive estimations of the existing biodiversity and its evolutionary history.

A distinct physiological role of MutY in mutation prevention in mycobacteria

Oxidative damage to DNA results in the occurrence of 7,8-dihydro-B-oxoguanine (8-oxoG) in the genome. In eubacteria, repair of such damage is initiated by two major base-excision repair enzymes, MutM and MutY. We generated a MutY-deficient strain of Mycobacterium smegmatis to investigate the role of this enzyme in DNA repair. The MutY deficiency in M. smegmatis did not result in either a noteworthy susceptibility to oxidative stress or an increase in the mutation rate. However, rifampicin resistant isolates of the MutY-deficient strain showed distinct mutations in the rifampicin-resistance-determining region of rpoB. Besides the expected C to A (or G to T) mutations, an increase in A to C (or T to G) mutations was also observed. Biochemical characterization of mycobacterial MutY (M. smegmatis and M. tuberculosis) revealed an expected excision of A opposite 8-oxoG in DNA. Additionally, excision of G and T opposite 8-oxoG was detected. MutY formed complexes with DNA containing 8-oxoG: A, 8-oxoG: G or 8-oxoG: T but not 8-oxoG : C pairs. Primer extension reactions in cell-free extracts of M. smegmatis suggested error-prone incorporation of nucleotides into the DNA. Based on these observations, we discuss the physiological role of MutY in specific mutation prevention in mycobacteria.

Synergistic effects of UdgB and Ung in mutation prevention and protection against commonly encountered DNA damaging agents in Mycobacterium smegmatis

The incorporation of dUMP during replication or the deamination of cytosine in DNA results in the occurrence of uracils in genomes. To maintain genomic integrity, uracil DNA glycosylases (UDGs) excise uracil from DNA and initiate the base-excision repair pathway. Here, we cloned, purified and biochemically characterized a family 5 UDG, UdgB, from Mycobacterium smegmatis to allow us to use it as a model organism to investigate the physiological significance of the novel enzyme. Studies with knockout strains showed that compared with the wild-type parent, the mutation rate of the udgB(-) strain was approximately twofold higher, whereas the mutation rate of a strain deficient in the family 1 UDG (ung(-)) was found to be similar to 8.4-fold higher. Interestingly, the mutation rate of the double-knockout (ung(-)ludgB(-)) strain was remarkably high, at similar to 19.6-fold. While CG to TA mutations predominated in the ung(-) and ung(-)/udgb(-) strains, AT to GC mutations were enhanced in the udgB(-) strain. The ung(-)/udgB(-) strain was notably more sensitive to acidified nitrite and hydrogen peroxide stresses compared with the single knockouts (ung(-) or udgB(-)). These observations reveal a synergistic effect of UdgB and Ung in DNA repair, and could have implications for the generation of attenuated strains of Mycobacterium tuberculosis.

Determination of Spontaneous Mutation Frequencies in Measles Virus under Nonselective Conditions

There is a paradox between the remarkable genetic stability of measles virus (MV) in the field and the high mutation rates implied by the frequency of the appearance of monoclonal antibody escape mutants generated when the virus is pressured to revert in vitro (S. J. Schrag, P. A. Rota, and W. J. Bellini, J. Virol. 73: 51-54, 1999). We established a highly sensitive assay to determine frequencies of various categories of mutations in large populations of wild-type and laboratory-adapted MVs using recombinant viruses containing an additional transcription unit (ATU) encoding enhanced green fluorescent protein (EGFP). Single and double mutations were made in the fluorophore of EGFP to ablate fluorescence. The frequencies of reversion mutants in the population were determined by measuring the appearance of fluorescence indicating a revertant virus. This allows mutation rates to be measured under nonselective conditions, as phenotypic reversion to fluorescence requires only either a single-or a double-nucleotide change and amino acid substitution, which does not affect the length of the nonessential reporter protein expressed from the ATU. Mutation rates in MV are the same for wild-type and laboratory-adapted viruses, and they are an order of magnitude lower than the previous measurement assessed under selective conditions. The actual mutation rate for MV is approximately 1.8 x 10(-6) per base per replication event. Copyright © 2013, American Society for Microbiology. All Rights Reserved.